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ABSTRACT 

We measure the angular correlation function of radio galaxies selected by the 843 MHz 
Sydney University Molonglo Sky Survey (SUMSS). Wc find that the characteristic im- 
print of large-scale structure is clearly detectable, and that the survey is very uniform. 
Through comparison with similar analyses for other wide-area radio surveys - the 1400 
MHz NRAO VLA Sky Survey (NVSS) and the 325 MHz Westerbork Northern Sky 
Survey (WENSS) - we are able to derive consistent angular clustering parameters, 
including a steep slope for the clustering function, w{0) oc 9^^'^. We revise upwards 
previous estimates of the NVSS clustering amplitude, and find no evidence for depen- 
dence of clustering properties on radio frequency. It is important to incorporate the full 
covariance matrix when fitting parameters to the measured correlation function. Once 
the redshift distribution for mJy radio galaxies has been determined, these projected 
clustering measurements will permit a robust description of large-scale structure at 
z ^ 0.8, the median redshift of the sources. 

Key words: large-scale structure of Universe - galaxies: active - surveys 



1 INTRODUCTION 

[ Recent deep, wide- area radio surveys have made it possible 
. to measure the patterns of angular clustering in the radio 
' sky. Radio selection provides a window on the galaxy distri- 
bution that differs significantly from optical wavebands. Ra- 
dio sources probe the high-redshift universe: the broadness 
of the radio luminosity function ensures that the detected 
distribution of objects has a median redshift z ~ 0.8, inde- 
pendently of fiux-density threshold (Condon 1989). The pro- 
jected clustering signal originates to a significant extent from 
redshift z and serves as a robust tracer of structure at that 
epoch. However, the wide redshift range spanned carries a 
significant disadvantage: the angular clustering is vastly di- 
luted by the superposition of unrelated redshift slices and 
hence the residual signal is faint, an order of magnitude 
smaller than the projected clustering of the faintest Sloan 
galaxies. Moreover, the apparent brightness of a radio galaxy 
provides no indication of its radial distance, whereas in opti- 
cal wavebands galaxies have a characteristic luminosity that 
may be used to infer rough redshift information. 

Radio selection is independent of Galactic and inter- 
galactic extinction, so that reliably complete catalogues can 
be obtained over large areas. Even so, measuring the faint 
imprint of clustering has only been possible in the latest gen- 
eration of surveys: Faint Images of the Radio Sky at Twenty 
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centimetres (FIRST, Becker, White & Helfand 1995), the 
Westerbork Northern Sky Survey (WENSS, Rengelink et 
al. 1997) and the NRAO VLA Sky Survey (NVSS, Condon 
et al. 1998). Each of these databases lists > 10^ objects 
over ~ 10, 000 deg^ to limiting flux densities S'i4oo mhz '~ 3 
mJy. The surveys differ in observing frequencies and angular 
resolution, thus inter-comparisons can reveal subtle instru- 
mental and selection effects in individual catalogues (Blake 
& WaU 2002b). 

A recent addition to this complementary suite of sur- 
veys is the 843 MHz Sydney University Molonglo Sky Sur- 
vey (SUMSS). This is the first deep wide-area radio survey 
mapping the southern skies; and it does so at a frequency 
intermediate between that of WENSS (325 MHz) and NVSS 
(1400 MHz). Large-scale structure in the WENSS, FIRST 
and NVSS surveys has been quantified by various authors: 
Cress et al. (1996), Magliocchetti et al. (1998), Rengelink 
& Rottgering (1999), Blake & Wall (2002a, 2002b) and 
Overzier et al. (2003). In this study we apply a first clus- 
tering analysis to the SUMSS radio survey and compare the 
results with those derived from NVSS and WENSS, thus 
testing whether clustering is a function of radio frequency. 

It is worth emphasizing that clustering studies serve as 
a rigorous test of the calibration and large-scale uniformity 
of a catalogue. In particular, the angular correlation function 
w{d) analyzes the structure present as a function of angular 
scale 9. Instrumental effects in a radio survey typically man- 
ifest themselves on particular characteristic scales, and are 



2 Chris Blake et al. 




«-270° 

Figure 1. A plot of SUMSS catalogue entries with integrated fluxes exceeding 20 mjy in the region analyzed for large-scale structure. 
The outer circumference represents declination 5 = —50° and the centre of the diagram is 5 = —90°, circles being plotted every 10 
degrees. Galactic latitudes 6 = 0°, ±5° and ±10° are marked on the plot as respectively solid, dashed and dotted linos. 



usually rendered transparent by such an analysis. For ex- 
ample, unexpected features in the clustering pattern within 
the NVSS led to corrections to the data analysis pipeline for 
that survey (Blake & Wall 2002b). 

To de-project the angular clustering amplitude and re- 
cover the true spatial clustering properties of the sample, we 
need to know the redshift distribution of the sources, N{z). 
The radio source population at mJy flux-density levels is a 
mixture of nearby star-forming galaxies and more distant 
Active Galactic Nuclei (AGN) (Condon 1989). For example, 
approximately 5 per cent of NVSS sources are detected in 
the 2dF Galaxy Redshift Survey (Sadler et al. 2002), and 
40 per cent of these are star-forming. The 95 per cent of 
NVSS galaxies not associated with 2dFGRS galaxies are 
overwhelmingly radio AGN with redshifts z > 0.2, which 
implies that ~ 98 per cent of NVSS radio sources are dis- 
tant AGN. As their optical counterparts are often extremely 
faint, A^(2:) has not yet been measured for radio AGN at mJy 
flux-density levels. 

However, constraints on N{z) exist from luminosity- 
function models, and utilizing these (with large extrapo- 
lations) we can estimate the present-day clustering length 
ro of radio galaxies (e.g. Magliocchetti et al. 1998, Blake 
& Wall 2002a, Overzier et al. 2003). The values inferred, 
in the range ro — 7 10 Mpc, correspond to a clus- 
tering strength intermediate between normal galaxies and 
rich clusters, consistent with the nature of radio AGN as 
optically-luminous elliptical galaxies inhabiting moderately 



rich environments. An additional aim of this study is to es- 
tablish the "best value" of the projected clustering ampli- 
tude of radio galaxies, to be employed in a more rigorous 
derivation of the clustering length when N{z) is known. 



2 THE SUMSS CATALOGUE 
2.1 Summary 

The Sydney University Molonglo Sky Survey (SUMSS) is 
imaging the southern {5 < —30°) radio sky at an observing 
frequency of 843 MHz. The survey uses the Molonglo Ob- 
servatory Synthesis Telescope (MOST), a 1.6-km long cylin- 
drical paraboloid reflector located near Canberra in Aus- 
tralia. The angular resolution of the SUMSS mosaics is 45" x 
(45"cosec|5|). The survey design and science goals are de- 
scribed by Bock et al. (1999); observations began in 1997 
and completion is scheduled for the end of 2003. The flrst 
instalment of the SUMSS source catalogue was presented by 
Mauch et al. (2003), who showed that it is complete to flux 
density 5843 mhz ~ 8 mJy for declinations 5 < —50° and 
18 mJy for 5 > —50°. SUMSS has comparable sensitivity 
and resolution to the 1400 MHz NRAO VLA Sky Survey 
(NVSS) in the north (5 > -40°). 

This clustering study uses Version 1.0 of the SUMSS 
source catalogue, dated 25 February 2003, available online 
at www.astrop.physics.usyd.edu.au/sumsscat. This ver- 
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sion of the catalogue covers approximately 3500 deg'^ of the 
southern sky and contains 107,765 entries. We restricted our 
analysis to the region S < —50°, which has a significantly 
fainter completeness limit than S > —50°, and we placed 
additional angular masks around the locations of the Mag- 
ellanic Clouds; (5° < a < 15°, -75° < 5 < -65°) and 
(70° < a < 90°, -75° < <5 < -65°). Raw SUMSS images 
contain artefacts, grating rings and radial spokes (Bock et 
al. 1999), typically located near very bright sources. SUMSS 
identifies these image artefacts using a reliable decision tree 
algorithm; however, to be conservative, we placed circular 
masks of diameter 1° around catalogued sources brighter 
than 1 Jy. After these masking procedures, 68,373 catalogue 
entries remained with flux densities 5843 mhz > 8 mjy. This 
version of the SUMSS catalogue excludes fields located at 
Galactic latitudes |6| < 10°, thus we can assume that the 
vast majority of the objects are extragalactic in origin. The 
sky distribution of the 20 mJy sample is plotted in Figure 
1. 

The overlap region of SUMSS and NVSS (-40° <S < 
—30°) permits some comparisons and cross-checks. The frac- 
tion of SUMSS sources missing from the NVSS catalogue is 
« 0.5 per cent, the majority of which are artefacts left be- 
hind by the decision tree (these generally appear close to 
very bright sources which are masked in this study as de- 
scribed above). Furthermore, cross-matched radio sources 
have a median spectral index Si, oc i/"" * (Mauch et al. 2003 
Figure 8), consistent with previous determinations at ~ 1 
GHz. 

2.2 Uniformity of the SUMSS catalogue 

Spurious fluctuations in source surface density within a sur- 
vey cause systematic offsets in clustering measurements, un- 
less these fluctuations can be accurately modelled. Large- 
scale density gradients with a magnitude of > 5 per cent 
will drown out the cosmological signal in the angular cor- 
relation function on scales 6>1°, as we discuss in Section 
4.2. In a wide-area survey undertaken with an interferometer 
such as the Very Large Array, large-scale gradients typically 
appear as a function of declination (Blake & Wall 2002a), 
because the projection of the interferometric baselines on the 
sky inevitably changes as a function of declination, whereas 
in principle all right ascensions can be observed at transit. 
Furthermore, the "snapshot" radio images used to construct 
a survey such as the NVSS have relatively sparse iii;-plane 
coverage and require considerable processing and cleaning 
to remove artefacts. In contrast, the uw-coverage of SUMSS 
images is inherently continuous owing to the nature of the 
MOST: deconvolution is practically unnecessary and image 
cleaning is trouble-free (Bock et al. 1999). 

Figure 2 plots the source surface density of the SUMSS 
catalogue as a function of declination for flux-density thresh- 
olds 8 mJy, 15 mJy and 30 mJy. There are no apparent 
spurious gradients as a function of declination, although the 
scatter of the data points exceeds that predicted by Poisson 
statistics (the error bars). This excess variance is produced 
by clustering, by multiple-component sources and (poten- 
tially) by spurious fluctuations in surface density. We placed 
a quantitative upper limit on the rms spurious density fluc- 
tuations by assuming that these were entirely responsible 
for the excess variance. We modelled the surface density at 
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Figure 2. SUMSS source surface density as a function of decli- 
nation. Sources are binned in declination bands for flux-density 
thresholds 8 mJy (solid circles), 15 mJy (open circles) and 30 mJy 
(triangles). The dashed lines indicate the average surface density 
at each threshold. The error in the number of sources A' in a band 
is plotted as the "Poisson error" \/lV. Masked areas are excluded 
from this analysis. 

declination 5 as a sum of two components, a{5) = ao + cTi{5), 
where (Jq is a constant-density "unperturbed" level and 
cTi{5) is a fluctuating offset (such that oT = 0). The vari- 
ance of the measured density over the declination bins is 
the sum of a "Poisson" contribution (due to ctq) and a "sys- 
tematic" contribution (due to ai). We derived the variance 
of cri by subtracting the contribution of Poisson statistics 
from the total variance; the fractional rms fluctuation is 
then /flue — -y/Var((Ti)/(To. The resulting upper limits on 
spurious density fluctuations in the SUMSS catalogue are 
respectively 2.2%, 2.0% and 2.5% at thresholds 8 mJy, 15 
mJy and 30 mJy. 

A systematic variation in surface density a creates 
an approximately constant offset in 'w{9) of magnitude e'^, 
where e = {a — a)/a is the surface overdensity (Blake 2002). 
Thus we can write: 

AM^)=?=:^ = (/fl.e)^ (1) 

Our measurements detailed in the previous paragraph there- 
fore placed an upper limit Aw{9) ^ (0.02)^ = 4 x lO"** on 
the effect of source density gradients in SUMSS. Thus our 
determination of w{9) in Section 3 is not significantly af- 
fected by any such gradients. 



3 ANGULAR CORRELATION FUNCTION 
FOR THE SUMSS CATALOGUE 

We quantified the projected clustering in the SUMSS cat- 
alogue using the angular correlation function, 'w{9). This 
two-point statistic is determined by counting the number of 
pairs of catalogue entries, DD{9), as a function of angular 
separation 9, and comparing the result to the number of 
pairs expected by random chance. The value of w{9) is the 
fractional enhancement of pair counts above random. We 
computed wip) using the Landy-Szalay estimator (Landy 
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& Szalay 1993). Using the Hamilton estimator (Hamilton 
1993) led to identical results. The error in each separation 
bin is plotted in the Figures as a "Poisson" error: 



Aw{e) = 



(2) 



This is known to be a very good approximation for the vari- 
ance of this estimator of the correlation function. 

When fitting models to the measured w{9), we used 
the full covariance matrix of the separation bins, Cw The 
covariance matrix for a given correlation function model 
was derived using the analytic approximation of Bernstein 
(1994, equation 39), who presents a specific treatment of the 
Landy-Szalay estimator; see also Eisenstein & Zaldarriaga 
(2001). Bernstein's approximations - A'^ ^ 1, wiO) <^ 1, 
wn <^ w{9) - are valid for the present surveys. We made an 
additional approximation by neglecting higher-order corre- 
lations (q = 0). Once the correlation matrix has been com- 
puted, the chi-squared statistic is given by 



X 



(3) 



where Aw is a vector containing the differences between the 
measured values of wiff) and the model. 

Apparent clustering in radio source catalogues is due to 
two effects: the cosmological clustering of separate galaxies, 
and the resolution of individual radio galaxies into multiple 
components which appear as separate catalogue entries. A 
w(ff) analysis neatly separates these two contributions. Ex- 
cess pairs with separations &<0.1° are almost exclusively 
produced by multiple components; cosmological effects pre- 
dominate on scales 0>O.1°. A good fitting formula for the 
voiff) data is a sum of two power-laws with sharply differing 
slopes, resulting in a clear break at S ~ 0.1° (see Blake & 
Wall 2002a and Overzier et al. 2003 for further justification 
of this model). 

These effects are, as anticipated, present in the SUMSS 
catalogue. Figure 3 plots measurements of the SUMSS w(S) 
at flux-density thresholds 8 mjy (which is the completeness 
limit) and 15 mJy. At 8 mJy, a sum of two power-laws is not 
in fact a good fit to the data points (for Q > 0.02°), with 
a minimum reduced chi-squared statistic x?cd = 2.16 (de- 
rived using the full covariance matrix; the probability that 
should exceed this value by chance is Q = 1.0 x 10""^). 
However, at 10 mJy the fit has improved to x?cd — 1-66 
(Q = 0.024) and the corresponding value for 20 mJy is 
Xrcd — 1-37 (Q — 0.11). The poor fit at 8 mJy is proba- 
bly due to unidentified systematic effects near the complete- 
ness limit of the survey, and we used the more conservative 
threshold of 10 mJy for subsequent analysis (see Section 
4.3). 

We investigated the dependence on flux density of the 
two power-laws used to model w(ff). As suggested by Figure 
3 and quantifled in Section 4.4, we found no change with flux 
density in the amplitude of the shallower power-law, which 
dominates on scales 0>O.1° and reflects cosmological clus- 
tering. This finding is in agreement with other wide-area ra- 
dio surveys in the flux density range 3 mJy < 5*1400 mhz < 50 
mJy. Note that this result is in sharp contrast to clustering 
measured in optical wavebands, where the clustering ampli- 
tude dechnes rapidly as galaxies become fainter. The expla- 
nation for this difference is simple: unlike in optical wave- 
bands, fainter radio galaxies are not (on average) more dis- 
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Figure 3. Measurement of the SUMSS angular correlation func- 
tion at flux-density thresholds 5843 MHz = 8 rnjy (solid circles) 
and 15 mJy (open circles). The best- fitting sum of two power-laws 
for the 8 mJy data is overplotted. This is not a good fit, and we 
preferred to use a more conservative threshold of 10 mJy for the 
subsequent analysis. 



tant. As shown in Figure 3, the multiple-components piece of 
w(ff) (at small separations) has an amplitude that decreases 
with flux density. Blake & Wall (2002a) demonstrated that 
the expected dependence of this amplitude is I/ct, where a 
is the source surface density. 

At very small separations d < 0.02° there is a fall-off in 
the value of w{6) with decreasing 9. This is not real, but is 
a signal-to-noise effect caused by the failure of the survey to 
resolve weak double sources with separations slightly greater 
than the beam-width. At very large separations 9 ~ 10°, the 
value of w{9) is consistent with zero. This is strong evidence 
for a high degree of uniformity in the survey: as detailed in 
Section 2.2, calibration gradients produce spurious offsets in 
the value of w{9). 



4 COMPARISON WITH OTHER WIDE-AREA 
RADIO SURVEYS 

In this Section we compare the SUMSS clustering results 
with equivalent analyses of the 1400 MHz NRAO VLA 
Sky Survey (NVSS, Condon et al. 1998) and the 325 MHz 
Westerbork Northern Sky Survey (WENSS, Rengelink et al. 
1997). These different samples are summarized in Table 1. 
The three surveys have comparable angular resolutions but 
different observing frequencies, potentially probing different 
selections of the radio source population. It is of interest to 
discover whether clustering is a function of radio frequency. 
Cosmological clustering of radio galaxies is independent of 
flux density in the range 3 mJy < 5*1400 mhz < 50 mJy (Sec- 
tion 4.4), thus the choice of comparison flux-density thresh- 
old for NVSS and WENSS is not important. 



4.1 NVSS 

The angular correlation function for NVSS has been mea- 
sured by Blake & Wall (2002a) and Overzier et al. (2003). 
The result of Blake & Wall for flux-density threshold 
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Table 1. Summary of radio source samples used in this investigation. 



Survey 


Frequency 


Angular resolution 


Declination 


Angular area 


Flux density 


Number of 


Surface density 




(MHz) 


(arcsec) 


range 


(deg2) 


(mJy) 


sources 


(deg-2) 


SUMSS 


843 


45 X (45 cosec|5|) 


-90° <5 < -50° 


2,516 


> 8 


68,373 


27.2 












> 10 


57,071 


22.7 












> 15 


40,843 


16.2 












> 30 


22,301 


8.9 


NVSS 


1400 


45 


-40° < (5 < 90° 


31,104 


> 10 


524,808 


16.9 


WENSS 


325 


54 X (54 cosec|(5|) 


30° < (5 < 75° 


8,409 


> 35 (int) 


107,937 


12.8 








30° < (5 < 60° 


6,579 


> 35 (peak) 


87,501 


13.3 




Angular separation 8 / deg 

Figure 4. Measurement of the NVSS angular correlation function 
at flux-density threshold 51400 MHz = 10 mJy, from Blake & Wall 
(2002a). The best- fitting sum of two power-laws is overplotted. 
This is not a good fit owing to a deficit of pairs with separations 
e Si 0.2°. 

51400 MHz ~ 10 mJy is plotted in Figure 4; this is the deepest 
flux density for which systematic surface density gradients 
are negligible. A sum of two power-laws is an excellent fit to 
the data points, except for separations 0.1° < 9 < 0.3°, for 
which there is a significant deficit of pairs. This results in 
the overall best fit being poor (xLd = 2.21, Q = 7.0 x 10"*). 

We made the following improvement to previous anal- 
yses of the NVSS w{6): we suggest that the data points at 
6 « 0.2° are suspect. The NVSS synthesized beam features 
very strong, broad diffraction spikes (sidelobes) offset by 
0.2° —> 0.3° from the main lobe, with amplitudes which are 
almost 30 per cent of the peak (Condon et al. 1998 Section 
4.1 and Figure 17). These diffraction features correspond to 
the shortest baselines used in the VLA observations, which 
are about 37 metres (~ 170 wavelengths). 

It is very plausible that the required deep cleaning of the 
NVSS images causes the small deficit of galaxy pairs with 
separations 9 ~ 0.2°. The consequent CLEAN bias will not 
remove 10 mJy sources near stronger sources, but it might 
reduce the flux density of the fainter source in each 0.2° pair 
by some fraction of the rms noise, which is (Tnoisc = 0.5 mJy, 
independent of the true source flux density (J. Condon, priv. 
comm.). The observed shortfall of 0.2° pairs is ~ 1 per cent 
(Figure 4). In comparison, the number of data pairs in the 
separation range 0.1° < 9 < 0.3° is increased by 1.6 per cent 
when the flux-density threshold is reduced by anoisc/5 from 



10.0 mJy to 9.9 mJy; we conclude that the shortfall of 0.2° 
pairs is probably an instrumental effect. 

The effect of the deficit of pairs with separations 9 ~ 
0.2° is to reduce spuriously the amplitude of the fitted 
clustering power-law, which must attempt to accommodate 
these low points. Fitting the clustering power-law to the 
separation range 9 = 0.3° — > 10° instead, a much better fit 
to the data was obtained (x?cd = 1.07, Q = 0.38) with a 
clustering amplitude ~ 50 per cent higher than previously 
recorded. We postpone further discussion of the best-fitting 
clustering parameters to Section 4.3. 

Overzier et al. (2003) describe an "unexplained bump" 
at separations 0.1° < 6* < 0.3° in the NVSS w{9) measured 
above a threshold of 200 mJy. This is very probably an ad- 
ditional consequence of the diffraction sidelobes. 

4.2 WENSS 

The first measurements of the WENSS angular correla- 
tion function were presented by Rengelink & Rottgering 
(1999); we performed our own investigation. We analyzed 
the "main" WENSS catalogue, which covers the region 
0° < a < 360°, 30° < S < 75°, and is available online at 
www.strw.leidenuniv.nl/~dpf/wenss. The catalogue lists 
merged multiple-component ("M") sources as well as their 
separate components ("C" sources); we omitted the "M" 
sources. In addition we masked areas within 5° of the Galac- 
tic plane, and inside a large box (15° x 15°) surrounding the 
bright radio source Cygnus A. 

The completeness limit of WENSS has been determined 
as 5325 MHz ~ 35 mJy (Rengelink & Rottgering 1999). Fig- 
ure 5 plots the WENSS catalogue surface density at this 
limit as a function of declination in two different ways: firstly 
for a 35 mJy threshold in (deconvolved) integrated fiux den- 
sity, and secondly for a 35 mjy/beam threshold in (con- 
volved) peak fiux density. For a point source the integrated 
and peak fiux densities are equal; the large WENSS beam 
ensures that these quantities do not usually differ signifi- 
cantly for catalogue entries. 

For the threshold in integrated fiux density in Figure 
5, there is a systematic gradient amounting to a ~ 10 per 
cent fluctuation between 5 = 30° and S — 75°. A possible 
explanation for this feature is that at higher declinations, 
the projection of the interferometric baselines on the sky is 
shorter, implying reduced surface brightness sensitivity and 
marginally reduced completeness as regards integrated flux. 
This surface density gradient creates a signiflcant spurious 
offset in w{9), drowning out the cosmological signal at sep- 
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o WENSS > 35 mJy/bean 
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Figure 5. WENSS source surface density as a function of decli- 
nation at integrated flux-density threshold 35 mJy (soHd circles) 
and peak flux-density threshold 35 mjy/beam (open circles). The 
dashed lines indicate the average surface densities at these thresh- 
olds. The error in the number of sources A'^ in a band is v^lV. 



arations 6>1° (see Figure 6). These systematic efltects are 
much less serious for the sample selected by peak flux den- 
sity. Thus, motivated by Figure 5, we created a more uniform 
WENSS source catalogue by imposing a peak flux-density 
threshold of 35 mJy/beam, and restricting our analysis to 
the decUnation band 30° < 5 < 60° 80 per cent of the 
available area). 

Our measurements of the WENSS wiO) for these two 
catalogue selections are displayed in Figure 6. For the peak 
flux-density threshold, a double power-law again provides a 
reasonable fit (x?od = 1-50, Q = 6.4 x 10~^). However, the 
influence of surface gradients is still not negligible and the 
angular correlation function at large separations {6 ~ 10°) is 
not consistent with zero. Therefore the best-fitting cluster- 
ing power-law spuriously acquires a significantly shallower 
slope than that determined for SUMSS and NVSS. We there- 
fore added an additional parameter - a constant offset - to 
the WENSS angular correlation function fitting formula. Al- 
though this extra parameter does not significantly improve 
the statistic of the best-fitting model, it is well-motivated 
by the existence of the surface density gradients. 



4.3 Comparison of clustering amplitudes 

We compared the angular correlation functions for SUMSS 
(> 10 mJy), NVSS (> 10 mJy) and WENSS (> 35 
mJy/beam) (Figure 7). We used the 10 mJy threshold for 
SUMSS (rather than the completeness limit of 8 mJy) be- 
cause 10 mJy is the deepest flux-density level for which a 
power-law model proved a good fit to the data. 

We fitted a power-law w{9) = AO'" to the cosmolog- 
ical clustering part of each correlation function, using the 
full covariance matrix. In the case of WENSS, we introduced 
a third variable parameter, a constant offset, to accommo- 
date the surface gradients. The angular ranges to which the 
functions were fitted were 6 > 0.2° (SUMSS, WENSS) and 
e > 0.3° (NVSS). For SUMSS and WENSS, 9 > 0.2° safely 
exceeds separations for which multiple component sources 
are important (6l<0.1°). For NVSS, 9 > 0.3° also avoids 
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Figure 6. Measurement of the WENSS angular correlation func- 
tion. The open circles plot the result for integrated flux-density 
threshold S'325MHz = 35 mJy over declination range 30° < S < 
75°. A spurious offset is evident for separations due to 

a source surface density gradient. This motivated us to select a 
more uniform catalogue by imposing a 35 mJy/beam threshold in 
peak flux density over the area bounded by 30° < S < 60°. The 
w{0) result for this second selection is plotted as the solid circles, 
together with the best-fitting double power-law model. 
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Figure 7. Overplot of the angular correlation functions 
SUMSS (> 10 mJy, solid circles), NVSS (> 10 mJy, open 
cles) and WENSS (> 35 mJy/beam, triangles). 



the influence of diffraction spikes at 6 ~ 0.2°. Figure 8 plots 
contours of constant in the space of the clustering param- 
eters {A, a). The results for the three independent surveys 
are consistent; there is no evidence that clustering is a func- 
tion of observing frequency. 

The best-fitting clustering parameters for the three 
surveys are tabulated in Table 2. They differ somewhat 
from those previously reported for the NVSS, which were 
A = (1.0 ± 0.2) X 10"^ a = 0.8 ± 0.1 (Blake & WaU 2002a, 
Overzier et al. 2003; note that 9 is measured in degrees for all 
quoted values of the amplitude A) . These differences are due 
to ignoring the low w{9) data at 6' ~ 0.2° in the NVSS mea- 
surement. The best-fitting slope is now somewhat higher, 
a ~ 1.1, agreeing with the value originally reported for the 
FIRST survey by Cress et al. (1996). In fact, it is not sur- 
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Table 2. Results of power-law fits to the cosmological clustering part of the angular correlation function, w{9) = A9~", for SUMSS, 
NVSS and WENSS. In the case of WENSS a constant offset is also fitted to the function. The best fit is obtained by minimizing the 
statistic, incorporating the full covariance matrix. The errors in the parameters are derived by varying each in turn from the best-fitting 
combination (keeping the others fixed) and determining the variation for which Ax^ = 1, the appropriate 1-sigma increment when 
varying one fitted parameter. 



Survey 


Flux density 


Amplitude A 


Slope a 


Constant offset 




Q 




(mjy) 


(xlO-3) 




(xl0~^) 






SUMSS 


> 10 


2.04 ±0.38 


1.24 ±0.16 




1.79 


0.03 


NVSS 


> 10 


1.49 ±0.15 


1.05 ±0.10 




0.89 


0.57 


WENSS 


> 35 (peak) 


1.01 ± 0.35 


1.22 ±0.33 


0.35 ±0.19 


1.37 


0.15 



a 



SUMSS > 10 mJy 
NVSS > 10 mJy 
WENSS > 35 mJy/beam. 




1 2 3 

Amplitude A / x lO"'' 

Figure 8. Measurement of the clustering parameters w{9) = 
AO-" for SUMSS (solid line), NVSS (dashed line) and WENSS 
(dotted line). Contours of constant shown in the space of 

the parameters (A, a), corresponding to Icr and 2cr confidence re- 
gions. As the plot is in the space of two varying parameters, the 
la and 2a contours are defined by increasing by respectively 
2.30 and 6.17 from its minimum. When generating this plot, the 
covariance matrix was held fixed at that corresponding to the 
best-fitting correlation function model. 



prising that the clustering slope should exceed the canonical 
value a — 0.8 of normal galaxies; it is known that a depends 
on galaxy type and takes on its steepest value for elliptical 
galaxies such as radio AGN (Norberg et al. 2002). The new 
value of the clustering amplitude, A ~ 1.6 x 10~^, also ex- 
ceeds the previous determination. 

Inclusion of the full covariance matrix in the fitting pro- 
cedure increases the errors on the clustering parameters by 
~ 50 per cent, compared to the approximation of neglecting 
correlations between separation bins. This result holds de- 
spite the fact that the variances of each bin are dominated 
by Poisson noise rather than by cosmic variance. 

The multiple-components piece of w{9), at 6<0.1°, 
varies between the three surveys (Figure 7). These differ- 
ences are explained by the 1/cr dependence of the amplitude 
of the small-angle w{6) on surface density a, together with 
the small variations in angular resolution between the sur- 
veys. 
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Figure 9. Angular correlation function amplitudes for SUMSS, 
NVSS and WENSS for various independent fiux-density ranges. 
The slope of w(0) was fixed at —1.1 in all cases. 



4.4 Projected clustering amplitude as a function 
of flux density 

We repeated the angular correlation function measurements 
for SUMSS, NVSS and WENSS for various independent 
flux-density ranges, and fitted a power-law w{9) = A6~^'^ 
to the cosmological clustering parts of each function (keep- 
ing the slope fixed to permit comparison of amplitudes). We 
fitted an additional constant offset in the case of the WENSS 
analysis. The amplitudes are displayed in Figure 9; results 
at different observing frequencies have been transformed to 
1400 MHz assuming a spectral index 5*^ oc 

Figure 9 demonstrates that the amplitude of projected 
clustering does not depend on either observing frequency 
or flux density (for 51400 mhz ^ 50 m Jy) . At brighter flux 
densities the signal-to-noise of the measurement is too low to 
allow flrm conclusions to be drawn. We note that Overzier et 
al. (2003) found evidence for an enhanced angular clustering 
amplitude in NVSS for a threshold 5*1400 mhz = 200 mJy, 
with a large error bar. 



5 SUMMARY 

The results of this investigation may be summarized as fol- 
lows: 

• The angular correlation function for the SUMSS radio 
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survey is measured for the first time and is found to have the 
same double power-law structure as previously identified in 
NVSS. 

• Our analysis reveals no evidence for large-scale system- 
atic gradients in SUMSS, such as might be imprinted by 
instrumental effects or data reduction. 

• The NVSS w{0) is re-examined, and data points on an- 
gular scales 0.1° < 6 < 0.3° are excluded from the fitting 
process due to the possible influence of imperfectly-cleaned 
bright sidelobes. 

• The WENSS w(6) is re-examined, and a systematic sur- 
face density gradient as a function of declination is identified 
for sources selected above a fixed threshold in integrated flux 
density. A more uniform source catalogue can be created by 
imposing a threshold in peak flux density. 

• The parameters of the clustering power-law, w{0) = 
AO-", derived for SUMSS, NVSS and WENSS, are found 
to be in agreement, having values a ~ l.l and A ~ 1.6 x 
10~^. Errors in the parameters for the individual surveys are 
mapped in Figure 8 and Table 2. 

• Inclusion of the full covariance matrix of the separation 
bins significantly increases the errors on the fitted clustering 
parameters, compared to the assumption that the bins are 
statistically independent. 

• We find no evidence for a dependence of projected clus- 
tering amplitude on either flux density or observing fre- 
quency. 

The inter-comparison of these three wide-area radio surveys 
has resulted in the reflnement of our knowledge of the pro- 
jected clustering of radio galaxies. Once the redshift distri- 
bution of these sources has been determined at mjy flux- 
density levels, their angular clustering amplitude may be 
used to make a robust measurement of the spatial cluster- 
ing of luminous elliptical galaxies at « fs 0.8. 
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